arXiv:1509.03402vl [cond-mat.mes-hall] 11 Sep 2015 


An exactly solvable model for a strongly spin-orbit-coupled nanowire quantum dot 


Rui Lian-Ao Wu,^’^ Xuedong Hu,"^ and J. Q. You^A 

^ Quantum Physics and Quantum Information Division, 

Beijing Computational Science Research Center, Beijing 100094, China 
^Department of Theoretical Physics and History of Science, 

University of the Basque Country UPV/EHU, 48008, Spain 
^IKERBASQUE, Basque Foundation for Science, 48011 Bilbao, Spain 
‘^Department of Physics, University at Buffalo, SUNY, Buffalo, New York 14260, USA 

(Dated: September 14, 2015) 

In the presence of spin-orbit coupling, quantum models for semiconductor materials are generally 
not exactly solvable. As a result, understanding of the strong spin-orbit coupling effects in these 
systems remains poor. Here we develop an analytical method to solve the one-dimensional hard-wall 
quantum dot problem in the presence of strong spin-orbit coupling and magnetic field, which allows 
us to obtain exact eigenenergies and eigenstates of a single electron. With the help of the exact 
solution, we demonstrate some unique effects from the strong spin-orbit coupling in a semiconductor 
quantum dot, in particular the anisotropy of the electron g-factor and its tunability. 


Introduction. —Semiconductor materials with strong 
spin-orbit coupling (SOC) have attracted increasing at¬ 
tention because of their importance in both fundamen¬ 
tal science and potential practical applications. No¬ 
table examples include topological insulator phase [1] 
in strongly spin-orbit coupled HgTe/CdTe quantum 
wells [2], and Majorana-fermion edge states in a semi¬ 
conductor nanowire [3-5] in proximity to a conventional 
s-wave superconductor. Strong SOC has also been used 
to electrically manipulate spin qubits confined in semi¬ 
conductor quantum dots (QDs) [6-10]. Furthermore, 
the Rashba SOC [11] is tunable by an external electric 
field [12], which provides a versatile platform for explor¬ 
ing physical effects in the strong SOC regime. In previous 
studies, however, SOC is usually treated perturbatively, 
which is clearly not ideal if it is strong. An exactly solv¬ 
able theoretical model would no doubt help clarify the 
physical picture of the strong SOC limit [13-17]. 

The solution of a quantum model is determined by 
the corresponding Schodinger equation with appropriate 
boundary conditions. Recently, an exact solution to the 
quantum Rabi model (a two-level system coupled to a 
harmonic oscillator) was obtained from a boundary con¬ 
dition based on the Z 2 symmetry of the model [18, 19]. 
Moreover, one can map the quantum Rabi model into 
the problem of one electron confined in a one-dimensional 
(ID) harmonic QD with SOC and magnetic field. Con¬ 
sidering that the harmonic and hard-wall confinement 
potentials are two of the best known exactly solvable 
problems for a particle without SOC, it is natural to ask 
whether an exact solution for the hard-wall potential in 
the presence of SOC can be found, even though they ap¬ 
pear to be quite different: one potential has boundaries 
at definitive positions, while the other does not. 

In this Letter, we show that the problem of one elec¬ 
tron confined in a ID hard-wall QD with strong spin- 
orbit interaction in a magnetic field is exactly solvable. 
We expand the Hilbert space with the eigenvectors of a 


harmonic oscillator, and derive a recursion relation for 
the expansion coefficients of the electron wave function. 
By imposing the hard-wall boundary condition, we de¬ 
rive a set of transcendental equations to obtain eigenen¬ 
ergies of the electron. We then use this solution to study 
strong SOC effects in a semiconductor QD. In particular, 
we show that the strong SOC modulates the electron g- 
factor in a semiconductor QD, giving rise to a periodic 
response to the magnetic-field direction. Furthermore, 
the g-factor can be tuned by adjusting either the direc¬ 
tion of the magnetic field or the size of the QD. 

Nanowire QD with a Zeeman field perpendicular to the 
spin-orbit field. —Consider a conduction electron that is 
confined in a thin semiconductor nanowire QD, and sub¬ 
ject to both an internal SOC field and an external mag¬ 
netic field. We first study the case where the magnetic 
field is perpendicular to the spin-orbit field and parallel 
to the quantum wire. The model Hamiltonian for this 
case is [20] 

H = +aa^p-\- -\-V{x), (1) 

2me 2 

where p = —ihdx is the canonical momentum, me is 
the effective electron mass, a is the Rashba SOC coef¬ 
ficient [11], A = geRsB is the Zeeman splitting, with 
ge and hb being the effective electron g-factor and the 
Bohr magneton, respectively. For the QD confinement 
potential V(x), we consider a ID square-well potential 
of infinite barrier height, i.e., V{x) =0 for jxj < L and 
V{x) = 00 for ja;j > L, where L is the size of the QD. 
The nominal material for the QD is InSb, and the Hamil¬ 
tonian parameters are explicitly listed in Table I. 

Hamiltonian (1) exhibits a Z 2 symmetry, for which the 
combined operator is a conserved quantity, where 
V is the parity operator: = 0. Thus a^V 

and the Hamiltonian H have common eigenfunctions 
'k = ['I'i(x),'I' 2 (a;)]"'", where 1111 ( 0 ;) and 1112 ( 0 :) are the 
two spin components of a given common eigenfunction. 
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TABLE I. Parameters of an InSb QD 


me/m^ 

Xso*’ (nm) 

9e 

B (Tesla) 

L (nm) 

0.0136 

50-200 

-50.6 

0.3 

50 


® m is the electron mass 
Xso = hl(mf^a) 


The eigenvalue equation of the operator a^V leads to 

T'i(x) =±^2(-x), (2) 

where ± is from a^V = ±1. In the meantime, the 
Schrodinger equation iJd' = E'i> can be expressed in 
terms of d'l and ^^2 as 

«'i(a^) + y^2(x) =0. (3) 

In fact, the Schrodinger equation iJt[/ = E'i/ in the 
dt 2 (a;)]"’" space corresponds to two equations. 
However, only one [e.g., Eq. (3)] is independent, while 
the other can be derived by combining Eqs. (2) and (3). 

The solution of the original single-electron problem can 
now be obtained by solving Eqs. (2) and (3). We first 
choose a complete set of basis states, with which the two 
components 'I'i(a;) and ^' 2 ( 2 ;) can be expanded. For the 
QD model considered here, we find that it is essential 
to choose the eigenfunctions of the harmonic oscillator, 
^—h ^rUeUjQX^'^ (t>n{x) = {n+^)hjjQ(j>n{x), as the basis 
states. These states can be expressed as [21] 

(jinix) = 1/ (ttcCq) 'H„{x/ xo)e-xjp [-x^/(2a;o)] , (4) 

where 'Hnix/xo) are the normalized Hermite polynomi¬ 
als, and xo = y^h/{meUJQ) is the characteristic length of 
the harmonic oscillator. With this basis, the matrix ele¬ 
ment of the spin-orbit interaction is non-vanishing only 
when An = ±1, so that the total Hamiltonian is a sparse 
matrix, and the iterative equations for the expansion co¬ 
efficients have an antomatic cut-off. Note that Xq (or ujq) 
is not a model parameter of Hamiltonian (1), so that the 
energy spectrum of the QD model should not depend on 
xo (or Wo), c.f. Fig. 1. 

Withont loss of generality, we can expand 'I'i(a;) as 

00 

'^l(x) = (5) 

n—0 

where are the expansion coefficients to be determined, 
and the introduction of imaginary unit i is to ensure that 
dn is real. According to Eq. (2), we can write 'I' 2 (a;) as 

00 

^' 2 ( 2 ;) = ^ ±(-i)"d„(/)„(x), (6) 

n—0 

where we have used the parity property of the eigen¬ 
function of the harmonic oscillator, i.e., <()„(— 2 ;) = 


2m, 


+ ap — E 



Spin-orbit length (nm) Spin-orbit length (nm) 


FIG. 1. Energy spectrum of the QD obtained by solving 
Eq. (8). (a) The ground-state energy of the QD calculated 
using harmonic-oscillator expansion basis states with differ¬ 
ent characteristic lengths: xo = 25, 50, and 75 nm. Inset: the 
hard-wall confining potential of the considered QD. (b) The 
lowest four energy levels of the QD. 


(—1)"^„(2;). Substituting 411 ( 2 ;) and 412 ( 2 :) in Eq. (3) 
with the expansions given in Eqs. (5) and (6), we obtain 
the following recursion relation for the expansion coeffi¬ 
cients [22]: 

\/(n — l)ndn -2 + r]V^dn-i + [2n -I- 1 — 4e ± (—1)"'2^] 

xdn + riy^8(n + l)d„+i = (n + l){n + 2)dn+2, (7) 

where e = E/(hwQ), ^ = A/(/iwo), and rj = xq/xso, with 
2:so = h/(meOi) being the spin-orbit length. We now let 
do = 1 and introduce a new real variable x = di/dp. Now 
all the expansion coefficients can be determined as dn = 
dn(s, x), and both 4li_2(2;) are obtained in terms of these 
coefficients. The energy spectrum of the Hamiltonian (1) 
can then be straightforwardly determined by imposing 
the hard-wall boundary condition: 4'i 2 (±T) = 0. In 
terms of the coefficients d„, the boundary condition can 
be expressed as [22] 

do "Ho ~ d2'H2 + d4'H4 — dg’Ho • • • = 0, 

diHi — do'Ho -|- d^T-L^ — d^'Hy • • • = 0, (8) 

where Hn = H„(T/xo). With two equations for only two 
variables e and Xi the energy spectrum can be exactly 
solved from Eq. (8). In addition, using the property of 
the Hermite polynomial [21], we can obtain the recursion 
relation for Hn, 

\/n + lH„+i = '/2(L/xo)'Hn - y/nUn-i- (9) 

Equations (7)-(9) are our central results. The energy 
spectrum of the considered Hamiltonian is determined by 
solving two transcendental equations instead of solving a 
matrix eigenvalue problem. As mentioned above, a key to 
our solution is that with the harmonic-oscillator basis, we 
are able to obtain the exact recursion relation in Eq. (7). 

Figure 1 shows the energy spectrum of the problem 
as a function of the spin-orbit length 2 : 30 , with panel 
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(a) showing the ground energy while panel (b) showing 
the lowest four eigenenergies. Recall that by choosing a 
complete set of harmonic-oscillator basis states, we have 
artificially introduced an additional parameter, the char¬ 
acteristic length (frequency) of the harmonic oscillator 
2^0 (wo) [see Eq. (4)]. Since xq (wq) is not a model pa¬ 
rameter of Hamiltonian (1), the energy spectrum of the 
system should be independent of this parameter. Indeed, 
our results in Fig. 1(a) verify this interesting point, that 
choosing different xo leads to the same eigenvalues for 
the considered Hamiltonian. 

Nanowire QD with a general magnetic field .—We now 
turn to the general case where the magnetic-field direc¬ 
tion is arbitrary in the x-z plane. The model Hamiltonian 
takes the form 

tP' a 

H = - \-aa^p-\ -(ct“ cos 9 + sind) + V(x), (10) 

znie 2 


exhibits the Z 2 symmetry. However, as we show below, 
this model is still exactly solvable. 

The two components 4>i(x) and 'I' 2 (a;) of the wave 
function 4' satisfy the Schrddinger equations 



A 

+ ap + — sin0 — E 
A 

— ap —— sin 9 — E 


4/1 -I- -^4'2COS0 = 0, 

4/2 + ^4'i cosd = 0. (11) 


We again choose the eigenfunctions of a ID harmonic 
oscillator as the basis states. In the absence of Z 2 sym¬ 
metry, we have to introduce two independent series of 
expansion coefficients for 4'i(a;) and 4'2(a;), 


^' 1 ( 2 :) = ^ j”e„(/)„(a;), 4'2(x) = ^ z”/„</>„ (x), (12) 

n=0 n=0 


where 9 is the azimuthal angle of the in-plane magnetic 
field. In contrast to Hamiltonian (1) for the special case 
of a perpendicular field, the Hamiltonian here no longer 

I 


where e„ and /„ are the coefficients to be determined. 
Substituting the two components 4^^ (i = 1,2) in 
Eq. (11), we obtain two sets of recursion relations 


\/n(n - l)e „_2 -I- V^r]e„-i + (S„,e -|- 2^sin9)e„ + (-l)”2^cos6»/„ -I- + l)r]en+i = -\/{n + l){n + 2)en+2, 

\/vpn^^fn -2 + V^rj fn-i + {^n,e + 2^sin6»)/„ -I- (-l)"2,fcos6>e„ -I- \/8{n + Vfr] f„+i = -sj (n + l)(n -|- 2)/„+2, 


(13) 


where = 2n-|-1 — 4e. While these recursion relations 

are much more complex than the one for the special case 
of Hamiltonian (1) with Z 2 symmetry, we can follow the 
same procedure to solve them. Specifically, we let Cq = 1, 
and introduce three new real variables, ei/eg = Xi, 
/o/eo = X 2 : and fi/co = xs- Now all the expansion 
coefficients can be determined as functions of these three 
parameters together with the rescaled energy e = E jhujo. 
Cn = e„(e,xi,X 2 ,X 3 ) and /„ = /n(e, Xi, X 2 , Xa)- Com¬ 
pared to the simpler Hamiltonian (1) with Z 2 symmetry, 
we now need to introduce three new real variables xi 1 X 2 > 
and X 3 i instead of one new real variable x- 

The energy spectrum of the system can again be deter¬ 
mined by the hard-wall boundary condition 4'i 2 (± T) = 
0 , which gives rise to 


coT-Lq — 62^.2 + CihLi — eerhie ■ ■ 

• = o. 


eiHi — epHs + epHs — 67^7 ■ • 

• = o. 


foPio — 12^.2 + ffU/L — fe'Hs ■ ■ 

• = o. 


flTdl — fsHs + /s'Hs — /7’H7 ■ 

■• = 0 . 

(14) 


These four equations are exact, and we have only four 
unknowns, £, xi: X 2 i and X3- Thus the four variables 
can be exactly determined, and the energy spectrum of 
the system obtained. 

The general solutions derived here converge to the spe¬ 


cial case solution when we take 9 = 0. Specifically, the 
two sets of recursion relations (13) are reduced to just 
one, in the form of recursion relation (7), and the bound¬ 
ary condition (14) is reduced to Eq. (8). Clearly, while 
Z 2 symmetry simplifies the calculations, it is not an es¬ 
sential ingredient in our solution here. Instead, it is the 
choice of the harmonic oscillator eigenstates as expansion 
basis that is crucial to our exact solutions of the model 
Hamiltonians (1) and (10). This is different from the 
quantum Rabi model, where the Z 2 symmetry plays an 
essential role to exactly solve the model [18, 19]. 

With the exact solution to the general Hamiltonian we 
can study a QD with strong SOC. One of the key proper¬ 
ties of a confined electron is its g-factor, which is strongly 
affected by SOC strength, and possibly by the QD con¬ 
finement. With our nominal example of InSb, its bulk 
g-factor is isotropic and takes the value ge = —50.6. Fig¬ 
ure 2 shows our calculated g-factor, which can be defined 
as geS = (£’2 — Ei)/{plbB) of a single confined electron 
with Hamiltonian (10). This relationship should be valid 
as long as Zeeman splitting is much smaller than the 
orbital excitation energy. Figure 2(a) shows the effec¬ 
tive g-factor as a function of the direction of the applied 
magnetic field, with a clear anisotropy as the field ro¬ 
tates [23-25] . It is a periodic function of 9 with a period 
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Dot size L (nm) 

FIG. 2. Effective electron g-factor of the QD Hamiltonian 
(10). Here the applied magnetic field is 0.3 T and the spin- 
orbit length is chosen as Xso = 50 nm. (a) The g-factor as 
a function of the applied magnetic-field direction, (b) The 
g-factor as a function of the QD size L. 


of TT. In particular, for 0 = 7r/2 or 37r/2, the operator cr^ 
is conserved by Hamiltonian (10), so that the SOC effect 
is straightforward to obtain, and the electron g-factor 
is equal to the bulk value = —50.6. For 0 = 0 or tt, 
the magnetic field is perpendicular to the spin-orbit field, 
so that is no long conserved, resulting in the largest 
modulation to the electron g-factor, with (/eff ~ —38.6. 
Figure 2(b) shows the g-factor dependence on the QD 
size L. The SOC strength in the QD is characterized by 
the ratio L/xso- When L increases, reducing the orbital 
energy splittings, the spin-orbit induced level mixing be¬ 
comes stronger, which is equivalent to an increase in the 
relative SOC strength. Consequently, increasing L (but 
still within the range when the orbital excitation energy 
is larger than the Zeeman splitting) enhances the mod¬ 
ulation to the electron g-factor. In short. Fig. 2 shows 
that we can manipulate the electron g-factor by varying 
the applied magnetic-field direction and the nanowire QD 
size. 

Summary .—In this study we have obtained the exact 
energy spectrum of a ID hard-wall QD in the presence of 
both a strong SOC and an applied magnetic field. The 
key ingredient of our solution is the selection of the bare 
harmonic-oscillator eigenfunctions as the basis states for 
the single-electron Hilbert space, which allows us to ex¬ 


press the system Hamiltonian as a sparse matrix and ob¬ 
tain a recursive relation for the wave function that can 
be solved. With the help of this solution, we are able to 
study strong SOC effects in in a ID semiconductor QD 
beyond the perturbation limit. In particular, we demon¬ 
strate a strong anisotropy in the electron g-factor, and 
the tuning of the electron g-factor via external param¬ 
eters, i.e., the QD size and the applied magnetic-field 
direction. 
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